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In this contribution we consider the issue of singularity resolution within loop 
quantum cosmology (LQC) for different homogeneous models. We present results of 
numerical evolutions of effective equations for both isotropic as well as anisotropic 
cosmologies, with and without spatial curvature. To address the issue of singularity 
resolution we examine the time evolution of geometrical and curvature invariants that 
yield information about the semiclassical spacetime geometry. We discuss generic 
behavior found for a variety of initial conditions. Finally, we show that the modi- 
fications which come from Loop Quantum Cosmology imply a non-chaotic effective 
behavior in the vacuum Bianchi IX model. 
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I. INTRODUCTION 

In general relativity (GR) the singularity theorems of Hawking, Penrose and Geroch tell 
us that, under reasonable assumptions, singularities are generic. A spacetime is said to be 
singular if it is not geodesically complete, which may happen when some geometrical curva- 
ture invariants diverge. The expectation is that, by quantizing the gravitational degrees of 
freedom, namely, with a complete theory unifying gravity and the quantum, the singularities 
shall be resolved. 

Loop quantization (as in Loop Quantum Gravity) of the homogeneous, isotropic and flat 
Friedman-Robertson- Walker (FRW) cosmology coupled to a massless scalar field </>, can be 
exactly solved pp. For that model it was shown that: 

• The matter density operator p has an absolute upper bound and the expansion 9 is 
also bounded. One can conclude that curvature scalars do not diverge. This is a signal 
that a singularity is not present. 

• All states undergo a bounce and with this, the big bang is replaced by a big bounce. 

• The GR dynamics is recovered as we go away from the Planck scale, this means that 
we are recovering the original theory that we quantize. 
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• Dynamics of semiclassical states are well captured by an effective theory that retains 
information about the loop quantum geometry. 

• Semiclassical states at late times must have come from semiclassical states before the 
bounce. That is, semiclassicality is preserved across the bounce j2]-FJ] 

• With all these, one can conclude that the singularities are resolved: the geodesies are 
inextendible, and are well defined on the other side of the would be big bang. 

The fact that the effective theory provides an accurate description of the dynamics at 
the Planck scale has been strongly used to explore the anisotropic models. The effective 
theory is obtained from the quantum Hamiltonian operator by taking expectation values 
on appropriately defined states. The thus obtained effective Hamiltonian then generates 
the dynamics on a classical phase space. The solutions to the effective theory were shown 
in [5] to accurately describe the evolution of the expectation value of the observables in 
the quantum theory when they are considered on semiclassical states. Those results were 
extended to open, closed and flat FRW models with and without cosmological constant 
(see [6] for a review). Loop quantum cosmology (LQC) has been extended to the simplest 
anisotropic cosmological models, namely Bianchi I, II and IX [7HH]. But in none of these 
cases has the quantum theory been solved, even numerically. Then, in order to study these 
models at the semiclassical level, one generally assumes that the effective theory reproduces 
the solutions to the quantum theory when semiclassical states are considered. This is our 
working hypothesis, which is well justified by the results in the isotropic cases. It would be 
interesting to know whether the evolution of the semiclassical states reproduce the solutions 
which we get from the effective theory. From this point of view, the study of the effective 
theory can be seen as the first step in this direction. 

The new issues to consider in the anisotropic models are: is the bounce generic? We 
now have anisotropy/Weyl curvature, how does it behave near the singularity/bounce? Can 
we have different kind of bounces, say, dominated by shear er 2 ? Are the geometric scalars 
such as the expansion 6, the shear a 2 and density p absolutely bounded? The goal of this 
contribution is to answer these questions using the effective theory for Bianchi I which has 
anisotropies, Bianchi II that has anisotropies and spatial curvature and Bianchi IX which 
has all the features of Bianchi I, II and is, furthermore, spatially compact. Even more, the 
Bianchi IX model has a non trivial classical limit, in the sense that, vacuum Bianchi IX is 
chaotic in the classical theory and behaves like Bianchi I with Bianchi II transitions as one 
approaches the singularity. In all the cases that we shall consider, the matter content shall 
consist of a massless scalar field playing the role of internal time. In the case of Bianchi IX 
we shall also consider the vacuum limit. 

II. PRELIMINARIES 

In this section we briefly review the quantization of some cosmological models which 
include k=0 and k=l FRW and Bianchi I, II and IX models by using loop quantum gravity 
methods. Let us consider the spacetime as M = X x R where S is a spatial 3-manifold 
which can be identified by the symmetry group of the chosen model and is endowed with 
a fiducial metric °q a b and associated fixed fiducial basis of 1-forms °oj l a and vectors °ef. If 
S is non-compact then we fix a fiducial cell, V, adapted to the fiducial triads with finite 
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volume V . We also define Li which is the length of the ith side of the cell along °ej and 

1 /3 

V = LiL 2 L 3 . We choose for compact S, Li = V with i = 1,2, 3. 

In general relativity, the gravitational phase space consists of pairs {A l al E" 1 ) on £ where 
A l a is a SU(2) connection and Ef is a densitized triad of weight 1. Since all of the models in 
which we are interested are homogeneous and, if we restrict ourselves to diagonal metrics, 
one can fix the gauge in such a way that A % a has 3 independent components, c\ and Ef has 
3 independent components, pi, 

K = Y°< and Et = P -^^-q°e1 (2.1) 

where pi in terms of the scale factors a« are \pi\ = LiLjdjCik (i ^ j ^ k). Note that in 
isotropic cases, each of phase space variables has only one independent component. Using 
(c l ,pi) for anisotropic models, the Poisson brackets can be expressed as {c\pj} = SnG^Sj 
and for isotropic models, the Poisson bracket is {c, p} = 87^7/3 where 7 is Barbero-Immirizi 
parameter. With this choice of variables and gauge fixing, the Gauss and diffeomorphism 
constraints are satisfied and the only constraint is the Hamiltonian constraint 



Cn — N 



v 



h ' h "' (^-(i + 7 2 X4)+K 



d d x, (2.2) 



with N the lapse function, TLmatter = pV an d &ab the curvature of spin connection T l a 
compatible with the triads. 

To construct the quantum kinematics, we have to select a set of elementary observables 
such that their associated operators are unambiguous. In loop quantum gravity they are the 
holonomies h e defined by the connection A\ along edges e and the fluxes of the densitized 
triad E? across surfaces. For our homogeneous models we choose holonomies and p^. To 
have the corresponding constraint operator, one needs to express it in terms of the chosen 
phase space functions h e and p{. The first term, e^EfE*/ as in loop quantum gravity, 
can be treated by using Thiemann's strategy [TO] . 

rpai rpbj i 

6ljk ^/W = ? 2^tW ^ ° U ^ h ^ {h * ' V]Tk) (2 ' 3) 

where is the holonomy along the edge parallel to z-th vector basis with length \x and 
V is the volume, which is equal to \Z\piP2pli\- Note that jj, is arbitrary. Now, to define an 



operator related to the first term of Eq.(2.2), we can use the right hand side of Eq.(2.3) and 
replace Poisson brackets with commutators. To find an operator related to the curvature 
F^ b , for isotropic models and Bianchi I, one can consider a square Dy in the i — j plane 
which is spanned by two of the fiducial triads (for the closed isotropic model since triads do 
not commute, to define this plane we use a triad and a right invariant vector °£f)> with each 
of its sides having length p!^ Therefore, F^ b is given by 

' /in . . — I 



, ( h u ■ - 1 ,\ ■ ■ 



(2.4) 



Since in loop quantum gravity, the area operator does not have a zero eigenvalue, one can 
take the limit of Eq.(2.4) to the point where the area is equal to the smallest eigenvalue of the 
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area operator, A 2 = Ay/Sn^l*, instead of zero. Then, fa[cn = A. We take fa\ = faLi where fa is 
a dimensionless parameter and, by previous considerations, is equal to fa = \y/\pi\/ y\PjPk\ 
(' ' .1 / /••)• 

For Bianchi II and IX, we cannot use this method because the resulting operator is not 
almost periodic, therefore we express the connection A % a in terms of holonomies and then 
use the standard definition of curvature F* b . The operators corresponding to the connection 
are given by [8] 

c, = ^ . (2.5) 
fa 

Note that using this quantization method for flat FRW and Bianchi I models, one has the 
same result as the direct quantization of curvature F^ b , but for closed FRW it leads to a 
different quantum theory which is more compatible with the isotropic limit of Bianchi IX. 
We call the first method of quantization curvature based quantization and the second one 
connection based quantization. In Bianchi II and Bianchi IX models the terms related to 
the curvatures, F^ b and Q* b , contain some negative powers of Pi which are not well defined 
operators. To solve this problem we use the same idea as Thiemann's strategy. 



w"-'* - - W^W ^'^'-' w " 2}) • (2 ' 6) 

where fa is the length of a curve, £ G (0, 1) and j G |N is for the representation. Therefore, 
for these three different operators we have three different curve lengths (/i, where /i 
and jl can be some arbitrary functions of pi, so for simplicity we can choose all of them to be 
equal to //. On the other hand we have another free parameter in the definition of negative 
powers of Pi where, for simplicity, we take j = 1/2. Since the largest negative power of Pi 
which appears in the constraint is —1/4, we will take I = 1/2 and obtain it directly from 



Eq.(2.6), and after that we express the other negative powers by them. The eigenvalues for 



the operator (pil -1 / 4 are given by 



MV, Pl ,p 2 ,p 3 ) = ^l[p] / \ (2.7) 



with 



h(V) = ^V + V c - tJ\V-V c \, and V c = 27r 7 A£j. (2.8) 

By using these results and choosing some factor ordering, we can construct the total 
constraint operator. Note that different choices of factor ordering will yield different opera- 
tors, but the main results will remain almost the same. By solving the constraint equation 
Ch ■ = 0, we can obtain the physical states and the physical Hilbert space "H p h y s- As a final 
step, one would need to identify the physical observables, that in our case would correspond 
to relational observables as functions of the internal time (j). 

To test singularity resolution we will study some geometric observables: expansion 9, shear 
a 2 , curvature scalars and also volume of the universe V and matter density p, as relational 
observables in terms of <fi. 

Since working with full quantum theories of the models is difficult and, as shown in [5] 
for some models, the behavior of the effective or semiclassical equations, which are 'classi- 
cal' equations with some quantum corrections, are good approximations to the numerical 
quantum evolutions even near the Planck scale, we will work with the effective equations. 
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III. EFFECTIVE THEORIES 

The purpose of this section is to present the effective theories obtained for the different 
models considered. It has four parts. In the first one, we consider the flat isotropic model 
and discuss its main features. In the second part we recall the k—1 FRW model for both 
quantizations available. The third part deals with the Bianchi I and II models. The Bianch 
IX model is considered in the last part. 



A. Isotropic Flat Model 

In the FRW model with k=0, the effective Hamiltonian is given by 

where p<f, is the momentum of the field, V is the volume and its conjugate variable. They 
are related to the c and p variables by the equations V = p 3//2 , (3 = cj ■ s fp and satisfy the 
Poisson bracket V} = AirG^ and {4>,p<p} = 1. It was shown [5] that the dynamics of 
semiclassical states are well captured by the effective Friedman equation 

H 2 = ^/i P 



Pcrit 



with H = V/3V the Hubble parameter and p = p^/2V 2 the matter density. The GR 
dynamics is recovered as we go away from the Planck scale, that is, for densities p < p crr t/10, 
with the critical density given by, p cr i t = 3/87rG7 2 A 2 ~ 0.41p p where p p is the Planck density. 
Thus, all trajectories undergo a bounce where the density reaches precisely its critical value 

Pcrit • 



B. Isotropic Closed Model 

Now, for the isotropic closed model, as we discussed in previous section, there are two 
different quantum theories depending on the two different methods of quantization of the 
curvature F* b . The Hamiltonians are [TTj, 

w£i = sJ^ sin2 ( A/3 - D > - sin2 D + & + 7 ^ 2 ] - § w ( 3 - 2 ) 

2 2 

= SvrlvA 2 ^ A/? " 2D Sin A/? + (1 + l2)E>2] ~ P ~2 * ° ' (3 - 3) 

where D = Xfi/V' 1 / 3 and •§ = (27T 2 ) 1 / 3 . From the Hamiltonians, the corresponding quantum 
corrected Friedmann equations are 



r2 87rG /„ .\A P-P2 



^^^(P-POII-^ 1 ) (3-5) 

Pcrit 
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where pi = p cri t[(l + 1 2 )D 2 - sin 2 /)], p x = p cvi a 2 D 2 and p 2 = p C ritD[(l + ^ 2 )D - 2 sin A/3]. 
Since for both effective theories there are some geometric observables which are not abso- 
lutely bounded, we go further and use more corrections which come from the inverse triad 
term in the full theory, to see if the unboundedness of those observables are generic, or 
whether they improve by adding more corrections. Therefore, the Hamiltonian constraints 
change to 

n ^ = [sin2(A/3 ~ D) - sin2 D + (1 + l2)D2] ~ y w (3 - 6) 

Wgi = W - 2D sin A/3 + (1 + ^D 2 ] - f * , (3.7) 

with 

A(V) = ^(V + V c -\V-V c \) = $y \I V < I < \ (3.8) 
is a correction term which comes from the operator e l ^E^E b J \/\q\. 



Bianchi I and II 



The effective Hamiltonian for Bianchi I and II can be written in a single expression, 
P1P2P3 



H 



BIl 



87rG7 2 A 2 
1 

+ 



sin piCi sm P2C2 + sin P2C2 sin p 3 c 3 + sin p 3 c 3 sin /X1C1 



8ttG 7 2 



a(p 2 P3) 3/2 . _ 

sm p\C\ 



0. 



The parameter a allows us to distinguish between Bianchi I (a = 0) and Bianchi II (a = 1). 
This Hamiltonian together with the Poisson Brackets {c l ,pj} = SirG^fS^ and {0, = 1 
gives the effective equations of motion. It has been shown that the density, the expansion 
and the shear are absolutely bounded for the case of Bianchi I [T2l [13] . The case of Bianch 
II is more involved, as we shall see later on. 



D. Bianchi IX 

In the previous cases, the effective Hamiltonian were found by choosing the lapse function 
as N = V. But now in Bianchi IX, we choose N = 1. This will allow us to include more 
inverse triad corrections. Then, the effective Hamiltonian is given by [2] 

V A A(V)h\V) ( . . . 



p\p\ sin p 3 c 3 + p\p\ sin p x c x + p\p\ sin p 2 c 2 



$A(y)h\V) 
' 47rCV r c 4 7 2 A 

h Q (V)V 2 2 



(P1P2) 4 + (PiPs) 4 + (P2P3) ' 
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IV. EFFECTIVE VACUUM BIANCHI IX 

In this section we study the vaccum Bianchi IX due to its non trivial classical behavior, 
in the sense that near to the classical singularity it presents a chaotic behavior described 
by the Belinski, Khalatnikov and Lifshitz (BKL) scenario. This case can be considered 
within our analysis because it is included in the Bianchi IX case with massless scalar field 



when we take the p^ = limit in Eq.(3.9). To study the BKL scenario within the effective 
theory we use the same idea given by Misner [151 IB], he., we study the potential term in 
the Hamiltonian constraint. Due to the new features introduced by the effective theory, to 
study the potential term is not enough. Therefore, we analyze also the density in order to 
determine how the BKL scenario changes. 



The Effective Potential 



It is helpful to use the potential term of the constraint to study the solutions. The classical 
potential which comes from the spin connection's curvature in the classical constraint, in 
terms of Misner variables is 



W = \e 



^e" 4/3 + - 4e _/3+ cosh V3f3_ + 2e~ w+ [cosh 2^/3- - 1]^ 



(4.1) 



where Q = — \ log V and the anisotropies /3± are defined via 



Oi 



-fi+03+W3/9_)/2 



a 2 = e -n+03 + -^-)/2 j a ^ 



(4.2) 



Since the Q dependence factorizes, one can obtain an anisotropy potential V(j3 + ,/3_) which 
exhibits exponential walls for large anisotropies. The universe can be seen as a particle 
moving in such a potential (W) that presents reflections at the walls. An infinite number 
of these reflections implies that the system behaves chaotically. When the volume becomes 
small, the quantum effects become important and one should work with the full quantum 
theory, but one can use the effective equations to have a qualitative view of what happens 



near the classical singularity. From the effective Hamiltonian Eq.(3.9), the modified potential 
can be derived as a function of Pi 



V 2 A(V)h 4 (V), , ., , 

w eS = y 4 - ( pI+pI+pI 



(P1P2) 4 + (piPs) 4 + (P2PsY 



2VV?J 



(4.3) 



For a simple case, when /3_ = and (3 + —> —00, the classical potential is W((3 + ,Q) ~ 
i e -4C-4/3+_ jf we rewr jt e the modified potential in terms of Misner variables, we can see 
that in this limit, the modified potential behaves as W e g ~ 2yQ 6 4/3+^ where the 
dependency of both classical and modified potential are the same. Thus, we have an infinite 
wall for the modified potential, too (See FigjTJ left). However, for small volumes, the modified 
potential can be negative at some points. 

One should note that the kinetic term is also modified. However, this modification does 
not change the qualitative behavior, because the kinetic term is always finite and different 
form zero for a finite time. 
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FIG. 1: Left, modified potential when /3_ = 0. Right, maximum allowed density vs. p\ and p2 
where p3 = p\. Both in Planck units. 

B. Density 

In the general it can be seen in Pig jl] (right), the maximum allowed density which 

arises from the modified Hamiltonian, has two distinct disconnected regions with positive 
values. If we impose the weak energy condition and start the evolution within one region, 
the universe cannot reach the other region. These two regions have different dynamics. To 
study the vacuum Bianchi IX, we start from large volumes which lie in region B of Figjl] 
(right) and, as we go to smaller volumes we cannot reach zero volume because 'crossing' to 
region A is not allowed. Therefore, there is a smallest reachable volume in region B and, 
since very large anisotropies are not allowed near this smallest volume, and the modified 
potential is not too large there, then we have, at most, finite oscillations before reaching the 
bounce. On the other hand, in the internal region A, the anisotropies are very large when 
some of the pi are very small, and then the volume of the universe cannot be large enough 
to start the evolution from there. 

V. RESULTS 

We will now compare the results of the effective theories for the isotropic FRW k=0 
and k—1, diagonal Bianchi I, II and IX. For all of them, the matter content consists of a 
massless scalar field satisfying the Klein-Gordon equation. A good starting point to compare 
the results is to answer the questions that we posed in the introduction, 

• Is the bounce generic? Yes. All solutions have a bounce. In other words, singularities 
are resolved. In the closed FRW and the Bianchi IX model, there is an infinite number 
of bounces and recollapses due to the compactness of the spatial manifold. 

• How does anisotropy/Weyl curvature behave near the bounce? For Bianchi I and II 
the shear and the curvature scalars, far from the bounce, are monotonic and approach 
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FIG. 2: Evolution of shear a 2 and Ricci scalar R for Bianchi II. The initial conditions at t = are: 
fiiCi = tt/2 and pi = 1000. 
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FIG. 3: Evolution of shear a 2 and Ricci scalar R for Bianchi IX. The initial conditions at t = 
are: /ZjCj = 7r/2, px = 90, p2 = 80 and = 100. 



their classical values, but when they reach the region near the bounce they behave 
differently. In Bianchi I, they present one maximum which occurs at the bounce. In 
Bianchi II, they exhibit a richer behavior, because now they can be zero at the bounce 
or near to it, see Fig. |2j and have more than one maximum (for the shear there are up 
to 4 maxima, shown in Fig. |2j and for the scalar curvature up to 2 maxima, see [PT]). 
In Bianchi IX, if we restrict the analysis to one of the infinite number of bounces, it 
can be shown that anisotropy and curvature behave as in the Bianchi I or II cases. Of 
currect research is whether there are new behaviors [18] . When we consider more than 
one bounce in the Bianchi IX case, there appear new interesting behaviors, shown in 
Fig. [| and Fig. |3| 

• Can we have different kind of bounces, say, dominated by shear a 2 ? Yes, but only 
in Bianchi II and IX. in Bianchi I the dynamical contribution from matter is always 
bigger than the one from the shear, even in the solution which reaches the maximal 
shear at the bounce [T7] . 

• Are geometric scalars 6, o 2 and p absolutely bounded? In the flat isotropic model all 
the solutions to the effective equations have a maximal density equal to the critical 
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time 

FIG. 4: Evolution of expansion 8 for Bianchi IX. 



density, p crit = 3/87rG7 2 A 2 , and a maximal expansion 



I 2 

max 



6nGp. 



crit 



3/(2 7 A)) 



when p = p cr it/2. For the FRW k = 1 model, every solution has its maximum density 
but in general the density is not absolutely bounded. In the effective theory which 
comes from connection based quantization, expansion can tend to infinity. For the 
other case, expansion has the same bound as in the flat FRW model. However, by 
adding some more corrections from inverse triad term, Eqs. (3.6, 3.7), one can show 



that actually in both effective theories the density and the expansion have finite values. 
For Bianchi I, in all the solutions p and 9 are upperly bounded by its values in the 
isotropic case and a 2 is bounded by <r max = 10.125/(37 2 A 2 ) [19]. For Bianchi II, 9, a 2 
and p are also bounded, but for larger values than the ones in Bianchi I, i.e., there are 
solutions where the matter density is larger than the critical density. With point-like 
and cigar-like classical singularities [T7], the density can achieve the maximal value 
(p ~ 0.54p p ) as a consequence of the shear being zero at the bounce and curvature 
different from zero. For Bianchi IX the behavior is the same as in the closed FRW; 
if the inverse triad corrections are not included, then the geometric scalars are not 
absolutely bounded. But if the inverse triad corrections are included then, on each 
solution, the geometric scalars are bounded but there is not an absolute bound for all 
the solutions H"8lfT9l. 



There are others results that are important to mention, 

• Classical and effective solutions are equal far away from the bounce. 

• Bianchi I, II and therefore the isotropic case k=0 are limiting cases of Bianchi IX, but 
they are not contained within Bianchi IX. While the isotropic FRW k=l is contained 
within Bianchi IX only if the inverse triad corrections are not included, when they are 
included then the k=l universe is a limiting case, like the k=0 universe. 

• A set of quantities that are very useful are the Kasner exponents (in classical Bianchi 
I, the scales factors are ai = t hi , where ki are the Kasner exponents), because they can 
be used to determine which kind of solution one obtains. The Kasner exponents tell 
us about the Bianchi I transitions |20j (if they exist) and particularly in Bianchi IX, 
they are used to study the BKL behavior in the vacuum case. 
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• Some important solutions are locally rotational symmetric. This means that, at each 
point the spacetime is invariant under rotation about a preferred spacelike axis, for 
example, if a 2 = a^. These solutions are such as the one with maximal density and 
the vacuum limit (where all the dynamical contribution come from the anisotropies) . 

VI. CONCLUSIONS 

One of the main issues that a quantum theory of gravity is expected to address is that of 
singularity resolution. Loop quantum cosmology has provided a complete description of the 
quantum dynamics in the case of isotropic cosmological models and singularity resolution 
has been shown to be generic. A pressing question is whether these results can be generalized 
to anisotropic models. In this case we lack a complete quantum theory, but one can rely 
on the existence of an effective description, capturing the main (loop) quantum geometric 
features. In this contribution we have described the main features of such effective solutions, 
for a variety of anisotropic cosmological models, with and without spatial curvature. As we 
have seen, singularities seem to be generically resolved as the time evolution of geometrical 
scalars is well behaved past the would-be classical singularity. The big bang is replaced by 
a big bounce. 

We have also studied the behavior of a modified potential for the vacuum Bianchi IX 
model when quantum effects become important. We showed that the potential wall does not 
disappear and we have, potentialy, chaotic behavior near the classical singularity. However, 
if the weak energy condition holds and, if we start from large volumes and evolve the 
equations into small volumes, there will be a lower bound for volume within region B (FigjTJ 
right), and one does not reach region A (connected to zero volume). Since there are no 
large anisotropies near the smallest allowed volume, the solutions will not exhibit chaotic 
behavior. 

With the study of these anisotropic models, a question that still arises is whether this 
bouncing non-singular behavior is generic for inhomogeneous configurations. That is, are 
we a step forward toward generic quantum singularity resolution? 
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